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Abstract. In this letter we present the discovery of a very light planetary companion to the star Ara (HD 160691). The 
planet orbits its host once every 9.5 days, and induces a sinusoidal radial velocity signal with a semi-amplitude of 4.1 ins -1 , 
the smallest Doppler amplitude detected so far. These values imply a mass of m,2 sini=14M^ (earth-masses). This detection 
represents the discovery of a planet with a mass slightly smaller than that of Uranus, the smallest "ice giant" in our Solar System. 
Whether this planet can be considered an ice giant or a super-earth planet is discussed in the context of the core-accretion and 
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migration models. 
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03 ' 1 ■ Introduction precision better than 1ms 1 , After only a few weeks of opera- 
tion, it has discovered a first "hot-jupiter" (Pepe et al. 2004) or- 

The discovery of giant planets around other solar-type stars „ , ri miKvnc t, u i if • • j- i 

1 b r Jr biting the K dwarf HD 330075. The level of precision in radial- 

has opened the way to a new era of planetary research. The , .. , ..Anno • t 

r J velocity measurements achieved with HARPS gives now, for 

new worlds present a wide variety of orbital characteristics c ... fl . ■ - c ,., lU j . 

r J the hrst time, the possibility of lowering significantly the detec- 

and minimum masses, and 9 years after the first announcement .. ,. ... <.v. » • ,., ,., , 

J tion limit to the few-earth-mass regime, provided that the sig- 

(Mayor & Oueloz 1995), some of their properties are still defy- , . , , , . „ .,, .. , , , 

v J ^ ' r r J nal induced by stellar oscillations can be reduced with the use 

ing the theories of planetary formation. The increasing number „ . . , , t , . N 

° r J ° of an appropriate observing strategy (Bouchy et al., in prep.), 

of known systems is, however, giving the possibility to explore T ... , „ . ,. „ , . , 

J ,"5 In this letter we present the discovery of a ^14-Mm short 

their properties from a statistical point of view (e.g. Santos et . , n c , s . , , . .... . . 

r r r v b period (P~9. 5 days) extra-solar planet orbiting the star /i Ara, 

al. 2001; Zucker & Mazeh 2002; Udry et al. 2003; Eggenberger . . iji tU u-^ju i j 

' J bb b a star that was already known to be orbited by a longer period 

et al. 2004), and the observational and theoretical approaches . . , . , , onn ^ ^ -4.1.4.1. 1 

' rr giant planet (Butler et al. 2001). Together with the very low 

are now starting to converge (e.g. Trilling et al. 2002; Alibert • . „ r /A/r » tU * 1 <->nn/i\ ^ 

6 & v & & > mass companion to 55 Cnc (McArthur et al. 2004), these are 

et al. 2004; Ida & Lin 2004a). , t F , . . , . , . '\ . ™ 

' the only two sub-neptunian planets discovered to date. They 

Recently, with rte installation of the ne^w HARPS spec- ^ suspected to be earth . Uke rocky p i anets , orbiting solar-type 

stars. 



trograph (Pepe et al. 2002) at the 3.6-m ESO telescope (La 
Silla, Chile) a significant quantitative advance has been pos- 
sible. This state of the art instrument is capable of attaining a 
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* Based on observations collected at La Silla Observatory, ESO, nitude southern G5V star in the constellation Ara, the Altar, and 

Chile, with the HARPS spectrograph, at the 3.6-m ESO telescope according to the Hipparcos catalog (ESA 1997), it has a paral- 

(programs 073.D-0578 and 072.C-0488). lax of 65.5±0.8 mas, which implies a distance from the Sun of 
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Fig. 1. HARPS radial- velocity measurements of n Ara as a function of 
time. The filled line represents the best fit to the data, obtained with 
the sum of a keplerian function and a linear trend, representing the 
effect of the long period companions to the system. The residuals of 
the fit, with an rms of only 0.9 m s" 1 , are shown in the lower panel. 



15.3 pc, and an absolute magnitude of M v =4.20. Its color index 
B - V is 0.694. 

From a HARPS spectrum with a S/N ratio of the order of 
^1000 (average of 275 individual spectra), we have derived the 
stellar parameters for \i Ara using a fully spectroscopic analy- 
sis (Santos et al. 2004). The resulting parameters (T e ff, log g, 
V t , [Fe/H])=(5813±40K, 4.25±0.07dex, 1.30±0.05kms-\ 
+0.32±0.05 dex), are in almost perfect agreement with the val- 
ues published in Santos et al. (2004), Bensby et al. (2004), 
and Laws et al. (2003). The surface gravity derived using the 
Hipparcos parallax and an effective temperature of 5800 K is 
4.25 dex (see e.g. Santos et al. 2004). 

Using the temperature, [Fe/H], absolute magnitude and 
bolometric correction (Flower 1996), we derived a stellar mass 
of 1.10±0.05Mq for /i Ara, from an interpolation of the theo- 
retical isochrones of Schaerer et al. (1993). This is in excellent 
agreement with the 1.08 and 1.14Mq derived by Butler et al. 
(2001) and Laws et al. (2003), respectively. Preliminary results 
from the asteroseismology analysis are also in excellent agree- 
ment with these values (Bazot et al., in prep.). 

From the width of the CORALIE Cross-Correlation 
Function (CCF) we have computed a projected rotational ve- 
locity of 2.4km s -1 for /i Ara (Santos et al. 2002). This value 
is in agreement with the low chromospheric activity level of 
the star, logi?^ 7{ =-5.034±0.006, obtained from the HARPS 
spectra. Similar values of —5.02 were obtained both from the 
CORALIE data (Santos et al. 2000) and by Henry et al. (1996) 
at different epochs. The inactivity of this star is further sup- 
ported by its low (and non-variable) X-ray luminosity (Marino 




Fig. 2. Phase-folded radial- velocity measurements of /i Ara after sub- 
traction of the linear trend shown in the upper panel of Fig. 1 . In both 
panels the error bars represent the rms around the weighted average of 
the individual measurements for a given night. 



2002), as well as by the lack of significant photometric varia- 
tion in the Hipparcos data (ESA 1997). 

From the observed value of log R' HK we can infer an age 
above ^2Gyr (Pace & Pasquini 2004) and a rotational pe- 
riod of ^31 days (Noyes et al. 1984). This age is compatible 
with the 4.5 Gyr obtained from an interpolation of theoretical 
isochrones (e.g. Laws et al. 2003), and with the upper value for 
the lithium abundance loge(Li)<0.86dex derived by Israelian 
et al. (2004) for this dwarf. 

3. Radial velocities 

In June 2004, /i Ara was intensively measured over 8 consec- 
utive nights with the HARPS spectrograph as part of an as- 
teroseismology program (Bouchy et al., in prep). During each 
night, we obtained more than 250 spectra of this star, from 
which we derived accurate radial velocities. The average radial 
velocity for each night was then computed from a weighted av- 
erage of each individual value, its precision being limited by 
the uncertainty in the wavelength calibration 1 . 

The main motivation of this program was to study the pos- 
sibility that the high metal content of the planet-host stars 
(e.g. Gonzalez 1998; Santos et al. 2001, 2004, and references 
therein) is due to the engulfment of metal rich planetary ma- 
terial into their convective envelopes. Although current stud- 
ies seem to favor that the observed "excess" metallicity re- 
flects a higher metal content of the cloud of gas and dust that 
gave origin to the star and planetary system, recent results have 
suggested that this matter may still be unsettled (e.g. Vauclair 
2004). The asteroseismology technique provides us with a good 
tool to possibly solve this problem. As shown by Bazot & 
Vauclair (2004), precise stellar oscillation measurements may 
be able to determine if there is some metallicity gradient in the 
stellar interior, that could be a hint of strong stellar "pollution" 



The nightly average of the HARPS radial velocities will be avail- 
able in electronic form at CDS 
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Table 1. Orbital elements of the fitted 9.5-days period orbit and main 
planetary properties. 
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* Equilibrium temperature computed with an albedo of 0.35. 

events. The results of the asteroseismology campaign will be 
presented in Bouchy et al. (in prep.) and Bazot et al. (in prep). 

A first analysis of the data revealed what could be a peri- 
odic variation with an amplitude of about 4ms _1 (see Figs. 1 
and 2). As part of the HARPS GTO program, this star was 
then closely followed from July 14th to August 19th 2004 (16 
radial-velocity measurements were obtained). Each night the 
radial velocity was measured from the average of about 15 con- 
secutive independent radial velocity estimates (computed from 
different spectra) taken during a period of ~20 minutes. This 
methodology makes it possible to average the radial-velocity 
variations due to stellar oscillations (Mayor et al. 2003) - see 
also Bouchy et al. (in prep.). As seen in Fig. 1, the measure- 
ments done during the first 8 nights (when the star was fol- 
lowed during the whole night) have a considerable lower rms 
around the best keplerian fit than the following measurements. 
This scatter results from the photon noise error (~20cms _1 ), 
the calibration uncertainty (~40cms _1 ), and from the stel- 
lar noise (^80cms _1 ) that is not completely averaged on the 
nights with only 15 radial velocity measurements (Bouchy et 
al., in prep.). 

/i Ara was previously announced to harbor a giant planet 
in a long period (^740days) orbit (Butler et al. 2001). This 
orbital solution has since been updated by Jones et al. (2002), 
who found that the residuals of the radial-velocity planetary fit 
followed a long term trend, due to the presence of a second 
body in the system. 

In Fig. 3 we plot the radial-velocity measurements of /i Ara 
obtained during the last 6 years using three different instru- 
ments (see figure caption), as well as the best 2-keplerian fit. 
The orbit of the ^740-day period planet (actually with a period 
of ~660 days) is confirmed. However, the orbital parameters of 
the second (longer period) companion are not well constrained; 
we find a strong degeneracy between the derived orbital period 
and the value of the orbital excentricity, making it possible to fit 
the data with the former parameter varying between ^3000 and 
10000 days. Although not precisely determined, the mass of 
this companion remains probably in the planet regime. Despite 
of the still unconstrained long period of this outer companion, 
some stability studies of the system has been discussed (e.g. 
Gozdziewski et al. 2003). 
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Fig. 3. Radial velocity measurements of fi Ara obtained during the past 
6 years with the CORALIE (dots) and HARPS spectrographs (open 
triangles), and by Jones et al. (2002) (open circles). The curve repre- 
sents the best 2-body keplerian fit to the data. In the lower panel we 
present the rms around the fit. For the longer period keplerian fit, the 
eccentricity was fixed to a value of 0.2. 

4. A 9.5-days period planet with 
14 Earth-masses 

In Figs. 1 and 2 we present the HARPS radial-velocity mea- 
surements of /i Ara as a function of time. In this figure, the 
curve represents the best fit to the data, obtained with the sum 
of a keplerian function and a linear trend. The derived slope of 
this trend is in agreement with the expected effect due to the 
longer period companions (see Fig. 3). 

The analysis of the radial velocity measurements reveals a 
variation with a period of 9.5 days, and a semi-amplitude of 
about 4ms" 1 . These values can be explained by the presence 
of a m2 sin i = 14 M ffl planet orbiting /i Ara in a circular orbit. 

The residuals around the best fit to the HARPS data are flat, 
with a rms of only of 0.9 m s _1 . This rms decreases to the cali- 
bration level (0.43 ms -1 ) for the first 8 nights, attesting the in- 
credible precision of this instrument. Despite the low amplitude 
of the radial velocity signal, the false alarm probability that it 
is due to random noise is lower than 1%, as derived through a 
Monte-Carlo simulation. 

From the stellar luminosity and effective temperature we 
can derive a radius of ^1.32 solar radii for fi Ara. Combined 
with the rotational period of 3 1 days (see Sect. 2), this implies 
a rotational velocity of the order of 2.2 km s _1 for /i Ara, close 
to the measured value v sini=2.4kms~ 1 . Supposing that the 
orbital plane is perpendicular to the stellar rotation axis, this 
means that the orbital inclination sin i is close to unity, and that 
the observed minimum mass for the planet is not very different 
from its real mass. 



4 



N.C. Santos et al.: The first 14 earth-mass exoplanet 



Using the HARPS spectra we have derived both an activ- 
ity index, based on Call H and K lines, and the bisector of 
the cross-correlation function from the individual spectra. No 
correlation is found between these quantities and the radial ve- 
locities within the measurement precision. Given the very low 
activity level of /LtAra and the inferred rotational period of 
~30 days, it is very unlikely that rotational modulation is ca- 
pable of producing the observed stable periodic radial-velocity 
variation. Furthermore, to have a rotational period of 9.5 days, 
this star would have to rotate at about 7kms~ 1 . Such a rota- 
tional velocity would imply a much younger age for \i Ara, not 
compatible with its low level of activity. 

The presence of a 14 M® planet around fi Ara thus remains 
the only credible explanation for the observed 9.5-days period 
radial-velocity variation. 

5. Discussion 

As current planetary formation models are still far from being 
able to account for all the amazing diversity observed amongst 
the exoplanets discovered thus far, we can only speculate on 
the true nature of the present object. 

First, given its location and the characteristics of the cen- 
tral star, it is unlikely that this object was in fact a much more 
massive giant planet which has lost a large fraction of its en- 
velope over its lifetime. This is supported by the fact that more 
massive planets exist orbiting much closer to stars with simi- 
lar characteristics and by calculations by Baraffe et al. (2004) 
and Lecavelier des Etangs et al. (2004) which show that only 
planets significanly less massive than Jupiter would evaporate 
at 0.09 AU. Except if outward migration has occurred, we con- 
clude that the mass of this object has always remained small. 

To understand the consequences of this, it is necessary to 
recall that in the current paradigm of giant planet formation, 
a core is formed first through the accretion of solid planetes- 
imals. Once this core reaches a critical mass (m cr ; t ), accre- 
tion of gas in a runaway fashion becomes possible and the 
mass of the planet increases rapidly (e.g. Ida & Lin 2004b). 
This therefore implies that the current object has never reached 
the critical mass, for otherwise the planet would have become 
much more massive. Furthermore, recent giant planet forma- 
tion models including disk evolution and migration (Alibert et 
al. 2004) have shown that these effects greatly shorten the for- 
mation time. Hence, it is unlikely that the planet has migrated 
over large distances before reaching its present location. It was 
thus probably formed inside the ice radius (^3.2 AU - Ida & 
Lin 2004a), and its composition should be dominated by rocky 
(telluric) material. We note that the high [Fe/H] of /i Ara makes 
this case possible (Ida & Lin 2004a). Curiously, with 14 M® 
and a=0.09 AU, this planet is near the borderline of the mass- 
period desert defined by Ida & Lin (2004b), where no planets 
are supposed to exist. 

The above considerations lead us towards the following 
scenario for the formation of the present planetary system. The 
more massive planet, with the present ^660 days period or- 
bit, begins to form first and migrates inwards while growing in 
mass. Towards the end of the lifetime of the disk, the smaller 
planet is formed inside the orbit of the larger one, probably at 



a distance not exceeding 3 AU. Thus, we expect this object to 
have a massive, essentially rocky core (as opposed to icy), sur- 
rounded by a gaseous envelope with ^5-10% of its mass. It 
therefore probably qualifies as a super-Earth and not as a failed 
ice-giant. 

The discovery of this extremely low-mass planet represents 
a new benchmark for planet surveys, and demonstrates the abil- 
ity of instruments like HARPS to detect telluric planets with 
just a few times the mass of the Earth. In the future these de- 
tections will give the possibility to study the low end of the 
planetary-mass distribution. This kind of planets may be rela- 
tively common, as according to recent simulations (Ida & Lin 
2004a), very low-mass planets may be more frequent than the 
previously found giant worlds. This is further supported by the 
recent detection of a first neptunian planet in a short period or- 
bit around 55 Cnc (McArthur et al. 2004) 2 . Such planets will 
be preferential targets for space missions like the photometric 
satellites COROT and Kepler. Furthermore, the discovery of 
such low mass planets around stars that have at least one more 
giant exoplanet, makes of these systems very interesting cases 
to understand the processes of planetary formation and evolu- 
tion. 
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